Magnetization patterns in ferromagnetic nano-elements as functions of complex 

variable. 
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Assumption of certain hierarchy of soft ferromagnet energy terms, realized in small enough flat 
nano-elements, allows to obtain explicit expressions for their magnetization distributions. By min- 
imising the energy terms sequentially, from most to the least important, magnetization distributions 
are expressed as solutions of Riemann-Hilbert boundary value problem for a function of complex 
variable. A number of free parameters, corresponding to positions of vortices and anti- vortices, still 
remain in the expression. These parameters can be found by computing and minimizing the total 
magnetic energy of the particle with no approximations. Thus, the presented approach is a factory of 
realistic Ritz functions for analytical (or numerical) micromagnetic calculations. These functions are 
so versatile, that they may even find applications on their own (e.g. for fitting magnetic microscopy 
images). Examples are given for multi- vortex magnetization distributions in circular cylinder, and 
for 2-dimensional domain walls in thin magnetic strips. 

PACS numbers: 75.60.Ch, 75.70.Kw, 85.70.Kh 
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Complex analysis is a natural language for expressing 
solutions of many physical problems^ like 2-dimensional 
electrostatics, planar flows of ideal incompressible fluid, 
2-d problems of theory of elasticity, heat flow and, more 
recently, many-body wave functionSiS, of electrons in- 
volved in the fractional quantum Hall effect. In all these 
cases, a certain sets of physical assumptions were identi- 
fied, permitting to map the complex calculus to a subset 
of physical problems in a particular domain. The goal of 
this work is to do the same for nano-magnetics, express- 
ing magnetic structures in flat cylindrical nano-elements 
via analytic functions of complex variable, employing the 
conformal mapping to account for element shapes. 

Consider a cylinder with an arbitrary face shape, 
shown in Fig. [T] together with Cartesian coordinate sys- 
tem r — {X, Y, Z}. Its face is denoted by V and its thick- 
ness is L. Supposing that the cylinder is made of soft fer- 
romagnetic material with saturation magnetization Ms 
and the exchange stiffness C, one can introduce the nor- 
malized magnetization vector rh{'r) = {'mx,TnY,'mz} = 
M{r)/Ms, \rh\ — 1 and, classically,— express its magnetic 
energy as 



e[m] 




'.dWi] ■ m I dV,(l) 



where Le = y/C/(jIoMj) is the exchange length (of 
the order of 20nm for typical soft magnets), V = 
{d/dX,d/dY,d/dZ}, hDlrh] is demagnetizing field, cre- 
ated by the magnetization distribution rn{r). Square 
brackets denote the functional dependence (on the whole 
function, as opposed to its value at a particular point). 
Both e and hjj are functionals of m. 

In the static case with no electric currents the mag- 
netic charge formalism^ can be used for calculating hu 
by, first, expressing it as a gradient of a scalar potential 
function ho = Vu(r), which is, in turn, a solution of the 




FIG. 1: Cylinder with Cartesian coordinate system axes. 



Poisson equation V'^u p with the requirement (due to 
the finite size of the particle) that both \r\u and |r|'^|Vu| 
are finite as |r| — >■ oo, p = V • to is the density of mag- 
netic charges (on the surface they are proportional to the 
normal component of M). 

Minimizing the energy functional ([1} one would recover 
both ground and metastable magnetization distributions 
in the particle. Unfortunately, the general analytical so- 
lution of the resulting system of coupled non-linear par- 
tial differential integral Euler equations is possible only 
in a very few cases (most notably for the magnetization 
states of small ellipsoidal particles)."* In general case ap- 
proximations are necessary. 

In the following, instead of minimizing all the energy 
terms simultaneously, the sequential minimization will be 
adopted. The energy terms are first sorted from the most 
important to the least important. We also start with the 
set of all possible functions. Then, in order of decreasing 
importance of terms, we sieve the current set of functions, 
keeping only the ones, minimizing (extremalizing, to be 
precise) the current term. The procedure is repeated with 
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the remaining functions and remaining terms. The ques- 
tion is then how far we can go and how many functions 
remain at the end ? The final functions, constructed in 
this way, are, obviously, influenced the most by the first 
considered (more important) terms. 

To reduce the dimensionality of the problem, let us re- 
strict consideration to the particles in the form of planar 
cylinders, by starting not from the set of all possible 3-d 
functions m(r), but from 2-d functions, assuming that 
the cylinder is thin enough that the magnetization dis- 
tribution is independent on Z, that is dfh[r)/dZ — 0. 
Then it is convenient to introduce the complex coordi- 
nate z = X + lY , I = V— 1- The set of functions can be 
sieved further by restricting the length of magnetization 
vector |m| — 1, deflning it via stereographic projection 

mx+imy = 2w(z, z)/(l + (2) 
mz = (1-Kz,z)|2)/(1 + Kz,z)|2), 

where line over the variable means complex conjugation 
{z = X — lY) and the newly introduced complex function 
of complex variable w(z,z) is not necessary holomorphic 
(that is, not necessarily differentiable) but just an ar- 
bitrary relationship between two complex numbers, as 
noted by its explicit dependence on z. 

To sort the energy terms, we first note that there are no 
magnetic monopoles and, thus, the total magnetic charge 
of the particle (consisting of volume charges, proportional 
to V • m, and surface charges, proportional to the normal 
component of m at the surface) is zero. As the particle 
is shrinking, the positive and negative magnetic charges 
move closer together and so their positive self-energy is 
more and more compensated by their negative interac- 
tion energy. Thus, in small enough particles exchange 
interaction is more important than magnetostatic. It is 
more important to minimize the exchange energy even at 
a cost of having some magnetic charges. 

Magnetic charges are also different. There are volume 
and surface ones. It can be deduced from a very general 
considerations that for small particles the surface effects 
dominate. Thus, the volume charges are the least impor- 
tant in our case. Noting that planar cylinders have two 
surfaces: face and side, let us also distinguish between 
two corresponding types of surface charges. Of these, 
the face charges are more important to reduce, because 
the face is bigger. 

Thus, let us assume the following order of sequential 
minimization: exchange energy, energy of face charges, 
energy of side charges, and, finally, the energy of volume 
magnetic charges. 

Introducing complex derivatives d/dz — [d/dX — 
id/dY)/2, d/dz = {d/dX + id/dY)/2 the exchange en- 
ergy density can be represented as 

E,^ ,9 8 { dw dw dw dw\ 

iVm,f= + . 3 

^ ^ (1 + WW)'' \oz oz az oz J 

Putting to zero the first variation of integral of this den- 



sity we get the following Euler equation 

d f dw\ 2w dw dw 
dz \dz J 1 + WW dz dz 

This equation is non-linear, it has several families of so- 
lutions (their complete set is still unknown). The first 
family was found by Belavin and Polyakov.— These solu- 
tions are called "solitons" and correspond to w = f{z) 
being an arbitrary analytic function (that is dw/dz = 0). 
This is obvious in complex notation, in the original work^ 
derivation was much more involved. There are much 
more particular solutions of this Euler equation,— but 
Belavin and Polyakov solitons are the only ones, having 
finite energy in infinite ferromagneti^ In restricted thin 
film geometry other solutions of this Euler equation start 
to become relevant. 

Such other family of solutions was discovered by 
David Gross^ they are called "merons" and expressed as 

w{z,'z) = f{z) I \J f{z)f{'z), which can be verified by di- 
rect substitution. Here f{z) is again an arbitrary analytic 
function. Merons have \w\ — 1 and, consequently from 
(12), TUz = everywhere. There are no magnetic surface 
charges on the cylinder faces (and so the magnetostatic 
energy is at absolute minimum,— zero) . The exchange en- 
ergy density of merons, unfortunately, has non-integrable 
singularities at zeros and poles of f{z), making the ex- 
change energy divergent for any non-trivial f{z). This 
divergence can be avoided at a cost of some face magnetic 
charges by joining solitons and merons continuously 

r fi^)/ei 1/(^)1 < ei 

«^(^.^) = <^ fiz)/^f{z)f{z) ei < |/(z)| < 62 , (5) 

[ /(z)/e2 |/(Z)| > 62 

for an arbitrary analytic function f{z) and two arbitrary 
real constants < ei < 62 < 00. Being composed of solu- 
tions of Euler equation, this function is locally extremal 
to the exchange energy functional, while the amount of 
face magnetic charges can be controlled by selection of 
the free constants Ci and 62. These constants^ allow to 
tighten meron arbitrarily close to singularities and zeroes 
of /(z) trading between the decrease of the magnetostatic 
energy of face charges and the increase of the exchange 
energy of the particle. 

The next energy term to sieve solutions through is 
the magnetostatic energy of side charges. Unlike the en- 
ergy of the face charges, this energy can be completely 
put to zero (its absolute minimum) without conflict with 
more important terms we have already minimized. The 
corresponding requirement-'^^ for f{z) is to find a func- 
tion f{z) analytical in the region T) in such a way that 

-R-S [/(C)'^(C)] = (no magnetization components normal 
to the side), where € C — dV is the boundary ofD, and 
rt(C) = nx{C,) + my{C,) is the complex normal to C. This 
is an instance of the well known linear Riemann-Hilbert 
problemji 

This problem and the resulting magnetization distri- 
butions were already considered in Refllll by reducing 



the Riemann-Hilbert problem to the problem of Hilbert 
Privalov on the miit disk, as described in ReflU. This pro- 
cedure, however, misses some of the physically relevant 
solutions. To solve the Riemann-Hilbert problem in this 
work the region T) is conformally mapped to the upper 
half plane i, such that Imi > 0. Then, the functions with 
non-essential singularities, satisfying condition of no nor- 
mal component to the boundary of the upper half plane 
(that is, real on the real axis) can be written, generally, 
as rational functions with real coefficients 



Em ,j 



(6) 



where gi and hi are arbitrary real numbers. 

In case of arbitrary shape of particle face, let us in- 
troduce the conformal mapping z = M(t), transforming 
the upper half plane t into the particle face z. Also note 
that roots of polynomial with real coefficients are either 
real or come in complex conjugate pairs. Since the roots 
of nominator and denominator correspond to zeros and 
poles of /(z), that is, to soliton "hats", covering meron 
singularities ([5]), it is convenient to express these polyno- 
mials directly in terms of their roots. Suppose that nom- 
inator (denominator) has (rip) pairs of complex roots 
at points a^, TTl {ci, q), such that Ima^ > (Imc^ > 0); 
as well as rrir (rir) real roots bj, Im bj = (dj, Imdj = 0). 
Then m = 2'mp + mr, n = 2np + nr and the function f{z) 
can be written parametrically as 



/ = M'{t) 
z = M{t). 



nr:''oK-^)(^"On;i^ofe-"0 
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This expression together with and ^ gives the fam- 
ily of Ritz functions for magnetization distributions, de- 
pending on a number of parameters a, &, c, d, ei, 62. 

Let us now see some examples. The case of magnetic 
disk is, probably, the simplest. The most general form of 
conformal transform of the upper half plane to the unit 
disk is 



M{t) = - 



e'°'{h + it) 
h — it ' 



(8) 



where h is an arbitrary real number and a € [0,27r). In 
case of m = 2, n = this solution coincides with the one, 
obtained in ReflTTI. and further, in a particular case of 
centered vortex, with an ansatz of Usov and Peschany.— 
For bigger m and n a particular multi- vortex magneti- 
zation distribution in a cylinder is shown in Fig. [2] with 
parameters, chosen for plausibility. 

The other interesting case is that of a 2-dimensional 
domain wall in thin strip. This problem had been treated 
numerically)^ but has resurfaced recently in connection 
with the idea of racetrack magnetic memory flS Noting 
that the conformal map from the upper half plane to the 
infinite strip < Imz < 1 is 



M{t) 
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FIG. 2: A multi- vortex magnetization distribution in circular 
cylinder, comparable to the numerical simulations of Reffl^. 
It is (0, (O, (0, dHI with h = A, a = -k/A, ei = 0.01, 62 = 
400, (aXcJ) = ({i,-4 + 4i,6 + 7i},{3},{2j,6 + 2i},{}). 
There are three vortices, two anti-vortices and two side-bound 
skyrmions, their centers, where the magnetization vector is 
vertical, are marked with dots. 




(9) 



FIG. 3: Some domain wall configurations in thin strip: 
a) vortex wall (comparable to the numerical simulation 
of Ref[3); b) anti- vortex wall; c) and d) two different 
arrangements of a pair of two vortex walls. As before, 
dots mark the singularities in the meron. It is ((2]), 
®, 0, ® with ei = 0.01, 62 = 400, and (a,6,c,d); 
a) ({/},{0},{},{5,-l/5}); b) ({}, {0, 5, -1/5}, {/}, {}); 
c) ({90i,j/90},{0},{},{360, 1/360, -20, -1/20}); d) 
({90i, i/90}, {0}, {}, {360, -1/360, -20, 1/20}). 



where additional parameters, related to a certain freedom 
to transform the upper half plane to itself, were omitted 
for clarity. Some of possible domain wall configurations 
are shown in Fig. [31 It is trivial to extend these ex- 
pressions to model the whole racetrack, containing many 
interacting domain walls of different types. 



4 



Concluding the step by step minimization procedure, 
we have sieved the set of all possible magnetization dis- 
tributions to obtain a much smaller set, parametrized 
by a finite number of scalar parameters. The values 
of these remaining parameters, corresponding to ground 
and metastable states of magnetization, can be found by 
minimizing the total energy ([T]). This problem is much 
simpler, compared to full solution of corresponding Eu- 
ler equations, and is tractable analytically in a number 
of interesting cases, e.g. like that of displaced magnetic 
vorte»i^ or of "C"-type low-symmetry magnetic states in 
circular cylinders.™ 

It is also possible to generalize these Ritz functions 
by relaxing the strict requirement of no normal magne- 
tization components introducing an additional parame- 
ter, corresponding to rescaling the particle. This works 
very well for treatment of quasi-uniform states in circular 
cylinders,^^ reproducing quite intricate measurements;^^ 
This generalization can be applied to strips (by rescaling 
their width) and should produce a very precise set of Ritz 
functions. 

There is no "real" unchangeable by evolution topolog- 
ical charge in finite particles, since the energy barrier, 
separating the states with different polynomial degrees 
of nominator or doniminator in ([T]), is finite (unlike the 
case of infinite ferromagnet, considered by Belavin and 
Polyakov). Nevertheless, since the exchange interaction 
is "the most important" in the particles of a few exchange 
lengths in size, the difference m — n has a pronounced ef- 
fect on the magnetic states and still remains useful for 
the purpose of classification even though there is no di- 
rect proportionality between the total energy and |to — n| 
as in Belavin-Polyakov case. 



While the exchange energy of magnetization distribu- 
tions (O, consisting both of soliton and meron parts, can 
be calculated analytically in general form by Greene and 
residue theorems, magnetostatic energy in ([T]) currently 
has to be calculated on a case by case basis. Finding 
generic expression for magnetostatic energy of magneti- 
zation distributions ([7]) is still an open problem. Ab- 
sence of such an expression and long-range character of 
dipolar interaction (resulting in substantial dependence 
of its energy on particle shape) makes it difficult to es- 
timate the relative importance of energy terms a priori 
for certain particle size. As with other approaches, in- 
volving trial functions, the rigorous statement is that the 
energy of these approximate solutions if always higher 
than that of exact ones. This allows to check the validity 
of obtained results a posteriori. As a rule of thumb, for 
particles with the dimensions of a few Le the magnetic 
states, their energies and various derivatives are usually 
in a good agreement to experiment. For larger particles 
results may vary. 

Functions ([7]) are so versatile that they might be use- 
ful on their own, beyond the Ritz method. For example, 
to enhance the resolution of magnetic force microscopy 
images of small particles by fitting the positions of topo- 
logical singularities, comparing calculated and measured 
forces on the tip. 

Magnetization dynamics can be considered 
separately2° on top of these distributions. 

The author would like to thank Andrei Bogatyrev from 
the Institute of Numerical Methods of Russian Academy 
of Sciences for his help with solution of Riemann-Hilbert 
problem, which have eventually led to the expression ^ . 
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